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Abstract. The coupled nonlinear volatility and option pricing model presented recently by 
Ivancevic is investigated, which generates a leverage effect, i.e., stock volatility is (negatively) 
correlated to stock returns, and can be regarded as a coupled nonlinear wave alternative of the 
Black-Scholes option pricing model. In this short report, we analytically propose the two-component 
financial rogue waves of the coupled nonlinear volatility and option pricing model without an 
embedded w— learning. Moreover, we exhibit their dynamical behaviors for chosen different pa- 
rameters. The two-component financial rogue wave solutions may be used to describe the possible 
physical mechanisms for the rogue wave phenomena and to further excite the possibility of relative 
researches and potential applications of rogue waves in the financial markets and other related fields. 
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I. INTRODUCTION 



The use of the term freak wave was first introduced to the scientific community by Draper [lj, |2J . Freak waves 
are also called rogue waves, giant waves, monster waves, killer waves, extreme waves, abnormal waves, huge 
waves, super waves, or gigantic waves, which are formed due to the interaction of the nonlinearity and dispersion 
in the wave propagation. Nowadays, there is no unique definite for the rogue waves, but, a wave called a rogue 
wave means that the wave height H (i.e., distance from trough to crest) or the crest height r\ c (i.e., distance from 
mean sea level to crest) exceeds a certain threshold related to the sea state Q. Generally speaking, the common 
criteria of rogue waves is that H > 2H S or n c > 1.5H S , where H s stands for the significant wave height Q. 
Rogue waves, as a special type of solitary waves, have drawn much attention in some fields of nonlinear science 

16 1 , Bose-Einstein condensations 0, [3] , atmospherics [l|3] > and even 



such as oceanics [2|-|7| , nonlinear optics 
finance The first-order rational solution of the self-focusing nonlinear Schrodingcr (NLS) equation was 



first found by Peregrine 2l| to describe the rogue wave phenomenon, which was known as Peregrine soliton 



(or Peregrine breather) and could be regarded as a limiting case of the Ma breather when the breathi ng p eriod 



approaches to infinity 22|, |24j or of the Akhmediev breather as the spatial period tends to infinity 
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More recently, multi-rogue wave solutions of the self-focusing NLS equation were also found on basis of the 
deformed Darboux transformation in Moreover, nonautonomous rogue wave solutions have also been 

found for the generalized NLS equations with variable coefficients in one-dimensional space [15L ll^ ] and in 
three-dimensional spaces [3] based on the similarity analysis idea. The term 'rogue (freak) waves' was called 
'rogons (frcakons)' if they reappeared virtually unaffected in size or shape shortly after their interactions [l5| . In 
nonlinear optics, optical rogue waves have been observed in two different experiments, based on two nonlinear 
optical models, i.e., the higher-order NLS equation Q and the NLS equation 0, respectively. Moreover, they 
have verified that optical rogue waves reveal an approach to exert maximal control over a nonlinear system with 
minimal effort and can be used to produce highly stable supercontinuum 0]. 

As is well known, the celebrated Black-Scholes (alias as Black-Scholes-Merton) option pricing model 
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was presented, based on the the geometric Brownian motion (i.e. the stochastic differential equation) dS = 
fj,Sdt + adW(t) satisfied by the stock (asset) price S and the Ito lemma Q, where C(S,t) is the value of 
European call option on the asset price S at time t, fi is the instantaneous mean return, a is the stock volatility, 
W is a Wiener process, and r is the risk-free interest rate. This model has drawn much attention and also 
opens the new area of study in financial mathematics or financial engineering. The Black-Scholes model can be 
widely used for valuing the pricing of European-style options, but it cannot be used for valuing other exotic types 
of options such as American Options or Asian Options as it cannot incorporate exercise features or any path 
dependencies 28-3(3]. 



More recently, based on the modern adaptive markets hypothesis due to Lo 3jJ, |32(, Elliott wave market 



theory [33l . [3_4j , and quantum neural computation approach 35( , Ivancevic presented a novel nonlinear option 



pricing model (called the Ivancevic option pricing model) [36| 



M(s,t) 
1 dt 
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(2) 



in order to satisfy efficient and behavioral markets, and their essential nonlinear complexity, where tp(S,t) de- 
notes the option-price wave function, the dispersion frequency coefficient a is the volatility (which can be cither 
a constant or stochastic process itself), the Landau coefficient j3 = 0(r, w) represents the adaptive market po- 
tential. Some periodic wave solutions of Eq. ([2]) have been obtained |36|. Recently, we first presented the two 
financial rogue wave solutions of the Ivancevic option pricing model (|2|) and illustrated their dynamics (20I , which 
may explain real financial crisis/storms (e.g., 1997 Asian financial crisis/storm and the current global financial 
crisis/storm). In addition, for the aim of including a controlled stochastic volatility into the adaptive Ivancevic 
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option pricing model ([2]), the adaptive, symmetrically coupled, volatility and option pricing model (with interest 
rate r and Hebbian learning rate c) has been presented, which represents a bidirectional spatio-temporal associa- 
tive memory 



3G, 



37| . In the following, we would like to study the financial rogue waves of the coupled nonlinear 



model. 



II. THE COUPLED NONLINEAR VOLATILITY AND OPTION PRICING MODEL 



By introducing a controlled stochastic volatility into the adaptive- wave model © , Ivancevic also presented the 



coupled nonlinear volatility and option pricing model in the form 
.da(S,t) ld 2 a(S,t) 
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(3b) 
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where 



dt 2 OS 2 

to generate a leverage effect, i.e., stock volatility is (negatively) correlated to stock returns 
Eqs. (|3"a)l and (|3b|) are called the volatility model and option pricing model, respectively, ip(S, t) denotes the option 
pricing wave function, playing the role of a nonlinear coefficient in the volatility model (|3aj) , a(S,t) denotes the 
volatility wave function, standing for a nonlinear coefficient in the option pricing model (|3b[) . Both processes thus 
evolve in a common self-organizing market heat potential, so they effectively represent an adaptively controlled 
Brownian behavior of a hypothetical financial market, (3{r,w) = r^2^ =1 vjjgj with the adaptation equation 



Eqs. (|3ei|) and f|3b[) without an embedded w— learning (i.e., for constant f3 = r-the 



Wj = -Wj + c\a\gj\il>\ 

interest rate), possessed the solitary wave solutions [36 

To the best of our knowledge, the analytical financial rogue waves of Eqs. (|3a| and (|3bj) were not reported 
before. In the following, we would like to investigate the two-component financial rogue wave solutions of Eqs. (|3aj) 
and (|3b[) without an embedded w— learning, which may play an important role to explain some occurred financial 
crisis/storms (e.g., 1997 Asian financial crisis/storm and the current global financial crisis/storm). 



III. FINANCIAL ROGUE WAVES (ROGONS) 

Here, based on the similarity analysis and the approach developed in Refs. [ill (l2 
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we show that 

the coupled nonlinear volatility and option pricing model (|3a[) and f|3b[) without an embedded w— learning also 
admits the financial multi-rogon (rogue wave) solutions, which may be used to describe the possible formation 
mechanisms for rogue wave phenomenon in financial markets. In the following, we only give the first two 
representative financial rogue wave (rogon) solutions of Eqs. Q3aj) and (|3bj) without an embedded w— learning . 

A. Financial one-rogue waves (rogons) 



The financial one-rogon solution of Eqs. pap and (|3b[) without an embedded w— learning for the coupled 
volatility wave function a(S,t) and option pricing wave function ?p(S,t) can be written as 

4(1 + ia 2 t) 



<ri(5,t) 



ab 
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(4b) 



^/2p(a 2 + b 2 ) V 1 + 2a 2 (5 - kt) 2 + aH 2 

in terms of the complex rational functions of the stock price S and time t, which involves five free parameters 
a, /3, a, b, and k to manage the different types of financial one-rogue wave propagations, whose intensity distri- 
butions \ai(S,t)\ 2 and \ipi(S,t)\ 2 are displayed in Figs. Q] and [2] for the chosen adaptive market potential {3 = 1, 
the scaling a = 1.5, amplitude parameters a = 2, b = 5, and the gauge k = 0, 0.5. Notice that time t in these two 
figures can be chosen to be negative since the solution is invariant under the translation transformation t — > t — to. 



FIG. 1. (color online). Financial rogue wave propagations for the coupled volatility wave function ai(S,t) and option 
pricing wave function tpi_(S,t) of the financial one-rogon solutions given by Eqs. (|4a[) and (|4b|) for parameters a — 1.5, /3 = 
1, a = 2, b = 5, and k — 0. (a) The intensity distribution \ai(S,t)\ 2 ; (b) the density distridution of \ai(S, t)\ 2 ; (c) The 
intensity distribution |cri (S 1 , £) | 2 for time t = (solid line), t = 0.4 (dashed line), t = 1.0 (dashed-dotted line); (d) The 
intensity distribution (S, t)\ 2 ; (e) the density distribution of \ipi(S,t)\ 2 ; (f) The intensity distribution \ipi (S, t)\ 2 for 
time t = (solid line), f = 0.4 (dashed line), t = 1.0 (dashed-dotted line). 
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FIG. 2. (color online). Financial rogue wave propagations for the coupled volatility wave function cri(S, i) and option 
pricing wave function ipi(S,t) of the financial one-rogon solutions given by Eqs. (|4a|) and (|4b[) for parameters a — 1.5, /3 = 
1, a = 2, b = 5, and fc = 0.5. (a) The intensity distribution \ai (S, t)\ 2 ; (b) the density distridution of \ai (S, t)\ 2 ; (c) The 
intensity distribution \ai(S, t)\ 2 for time t — (solid line), t — 0.4 (dashed line), t = 1.0 (dashed-dotted line); (d) The 



intensity distribution \ipi (S, t)\ 2 ; (e) the density distribution of \ipi(S,t)\ 2 ; 
time t = (solid line), t = 0.4 (dashed line), t = 1.0 (dashed-dotted line). 



(f) The intensity distribution \ip\{S,t)\ for 
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B. Financial two-rogue waves (rogons) 



Similarly, based on the symmetry analysis, we can obtain the financial two-rogon solutions of Eqs. pap and 
(|3b|) without an embedded w— learning for the coupled volatility wave function a(S, t) and option pricing wave 
function ip(S,t) in the form 

, P 2 (S,t) -(1/2) i a 2 t Q 2 (S,tf 



<J 2 (S,t) 

Ms,t) 



a a 



VWa 2 + b 2 ) 
ab 
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(5b) 



with these above-mentioned functions P 2 (S,t), Q 2 (S,t) and H 2 (S,t) being of polynomials in the stock price S 
and time t in the form 



P2(S,t) 



.£L ( S_ W )4_39L( S _ w) v 
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which contain five free parameters a, /3, a, &, and k to manage the different types of financial two-rogue wave 
propagations whose intensity distributions \a 2 (S,t)\ 2 and \ip 2 (S,t)\ 2 are depicted in Figs. |3] and |4] for the chosen 
adaptive market potential ft = 1, the scaling a = 1.5, amplitude parameters a = 2, b = 5 and the gauge k — 0, 0.5. 




FIG. 3. (color online). Financial rogue wave propagations for the coupled volatility wave function G2{S,i) and option 
pricing wave function if)2(S, t) of the financial two-rogon solutions given by Eqs. ()5a|) and (|5b|) for parameters a — 1.5, j3 — 
1, a = 2, b = 5, and k — 0. (a) The intensity distribution o^S, t)\ 2 ; (b) the density distridution of \a2(S, t)\ 2 ; (c) The 
intensity distribution \<T2(S, t)\ 2 for time t = (solid line), t = 0.4 (dashed line), t = 1.2 (dashed-dotted line); (d) The 
intensity distribution \ip2 (S,t)\ 2 ; (e) the density distribution of \ip2(S, t)\ 2 ; (f) The intensity distribution \%p2{S, t)\ 2 for 
time t = (solid line), f = 0.4 (dashed line), £ = 1.2 (dashed-dotted line). 
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FIG. 4. (color online). Financial rogue wave propagations for the coupled volatility wave function 02 (.?,£) and option 
pricing wave function ip^iS, t) of the financial two-rogon solutions given by Eqs. (|5a[l and (|5b[l for parameters a — 1.5, /3 = 
1, a = 2, b = 5, and k = 0.5. (a) The intensity distribution \tJ2(S, t)\ 2 ; (b) the density distridution of \o~2(S, t)\ 2 ; (c) The 
intensity distribution \(T2(S, t)\ 2 for time t — (solid line), t — 0.8 (dashed line), t = 1.5 (dashed-dotted line); (d) The 
intensity distribution \ip2(S, t)\ 2 ; (e) the density distribution of ^(S, t)| 2 ; (f) The intensity distribution \ip2(S, t)\ 2 for 
time f = (solid line), t = 0.8 (dashed line), t = 1.5 (dashed-dotted line). 

IV. CONCLUSIONS 

In conclusion, based on the symmetry analysis, we have investigated the coupled nonlinear volatility and option 
pricing model pap and f|3b[) without an embedded w— learning such that we present its analytical financial one- 
and two-rogon solutions with some free parameters. Moreover, we also illustrate their dynamical behaviors for 
chosen different parameters a, (3, k, a, and b (see Figs. QUI . 

Our results may play an important role to explain some real financial crisis/storms (e.g., 1997 Asian financial 
crisis/storm and the current global financial crisis/storm). Moreover, these results may further excite the pos- 
sibility of relative researches and potential applications for the financial rogue-wave phenomena in the financial 
markets and other related fields of science. 
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